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Abstract 

We calculate light-cone distribution amplitudes for non-relativistic bound states, including 
radiative corrections from relativistic gluon exchange to first order in the strong coupling 
constant. We distinguish between bound states of quarks with equal (or similar) mass, 
mi ~ m,2, an d between bound states where the quark masses are hierarchical, mi 3> 7712- 
For both cases we calculate the distribution amplitudes at the non-relativistic scale and 
discuss the renormalization-group evolution for the leading-twist and 2-particle distribu- 
tions. Our results apply to hard exclusive reactions with non-relativistic bound states in 
the QCD factorization approach like, for instance, B c — > rj c £u or e + e~ — > J/tpr] c . They also 
serve as a toy model for light-cone distribution amplitudes of light mesons or heavy B and 
D mesons, for which certain model-independent properties can be derived. In particular, 
we calculate the anomalous dimension for the B meson distribution amplitude 4>b(uj) in 
the Wandzura-Wilczek approximation and derive the according solution of the evolution 
equation at leading logarithmic accuracy. 



1 Also at: Tcchnische Universitat Miinchen, Physik Department, 85747 Garching, Germany. 
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1 Introduction 



Exclusive hadron reactions with large momentum transfer involve strong interaction dynamics at very 
different momentum scales. In cases where the hard-scattering process is dominated by light-like dis- 
tances, the long-distance hadronic information is given in terms of so-called light-cone distribution 
amplitudes (LCDAs) which are defined from hadron-to-vacuum matrix elements of non-local oper- 
ators with quark and gluon field operators separated along the light-cone [1,2] and [3,4]. LCDAs 
appear in the so-called pQCD approach to hard exclusive reactions [5-7], in the QCD factorization 
approach to heavy-to-light transitions [8], in soft-collinear effective theory [9,10], as well as in the 
light-cone sum rule approach to exclusive decay amplitudes [11-13] (for a recent review, see [14]). 

Representing universal hadronic properties, LCDAs can either be extracted from experimental 
data, or they have to be constrained by non-perturbative methods. The most extensively studied 
and probably best understood case is the leading-twist pion LCDA, for which phenomenological 
constraints [15-17] from the ir — 7 transition form factor [18], as well as estimates for the lowest 
moments from QCD sum rules [2, 19, 20] and lattice QCD [21, 22] exist. On the other hand, our 
knowledge on LCDAs for heavy B mesons [3,23,24], and even more so for heavy quarkonia [25,26], 
had been relatively poor until recently. 

Although LCDAs, in general, are not calculable in QCD perturbation theory, their evolution with 
the factorization scale (which is set by the momentum transfer of the hard process) can be calculated 
and is well understood, both, for light mesons [5] and for heavy mesons [23]. The situation becomes 
somewhat simpler, if the hadron under consideration can be approximated as a non-relativistic bound 
state of two sufficiently heavy quarks. In this case we expect exclusive matrix elements - like tran- 
sition form factors [27] and, in particular, the LCDAs - to be calculable perturbatively, since the 
quark masses provide an intrinsic physical infrared regulator. 

In this article, we are going to calculate the LCDAs for non-relativistic meson bound states 
including relativistic QCD corrections to first order in the strong coupling constant at the non- 
relativistic matching scale which is set by the mass of the lighter quark in the hadron. We discuss 
twist- 2 and twist-3 LCDAs for 2-particle Fock states with approximately equal quark masses (for 
instance an rj c meson), as well as 2-particle and 3-particle LCDAs for heavy mesons (like the B c ), 
where one of the quark masses is considered to be much larger than the second one (m& 3> m c ). Our 
results can also be viewed as a toy model for possible parameterizations of LCDAs for relativistic 
bound states, like the pion, kaon or B q meson at a low input scale, which may be evolved to the 
appropriate higher scales using the standard renormalization group equations in QCD (or HQET, 
respectively). 

Our paper is organized as follows. In the following section we give a short introduction to the 
non-relativistic approximation and collect the definitions and properties of LCDAs for light and 
heavy mesons. The main result of our paper, the corrections from relativistic gluon exchange, are 
presented in Section[3l Here we also derive model-independent results for the B meson LCDAs and 
namely the cut-off dependence of positive moments and the anomalous dimension kernels, and 
investigate the impact of the 3-particle LCDAs to the Wandzura-Wilczek approximation beyond tree- 
level. We discuss the effect of QCD evolution above the non-relativistic matching scale in Section [H 
including a new result for the B meson LCDA cf>g, before we conclude. Some technical details of the 
calculation are collected in an appendix. Some of our results have already appeared in a proceedings 
article [28]. 

2 Light-cone distribution amplitudes and the non-relativistic limit 

2.1 Non-relativistic approximation 

The wave function for a non-relativistic (NR) bound state of a quark and an antiquark with respective 
masses mi and mi can be obtained from the resummation of NR (potential) gluon exchange as 
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Figure 1: Rcsummation of potential gluons into a non-relativistic Coulomb wave- function. 



sketched in Figure[TJ The solution of the corresponding Schrodinger equation with Coulomb potential 
yields 



K 2 + \p\ 2 



2 ' 



(1) 



where k = m T a s Cp and m r = m\m 2 j{m\ +rri2) is the reduced mass. The normalization of the wave 
function gives the (non-relativistic) meson decay constant 



/nr. 



,3/2 



7T 



(mi + ?7T,2) 1 / 2 



(2) 



For more details and references to the original literature, see e.g. [29] (also [27]). 

In this approximation, the B c meson is entirely dominated by the 2-particle Fock state built 
from a bottom quark with mass m± = M = m^ and a charm antiquark with mass m 2 = m = m c . 
Consequently to first approximation in the NR expansion, the B c meson consists of a quark with 
momentum Mvu and an antiquark with momentum mv^, where Vn is the four- velocity of the B c 
meson (v 2 = 1). The spinor degrees of freedom for the B c meson are represented by the Dirac 
projector ^(1 + ^75. Similarly, a pseudoscalar rj c meson is interpreted as a cc bound state where 
both constituents have approximately equal momenta mv^. 

The non-relativistic approximation can also serve as a toy model for bound states of light (rel- 
ativistic) quarks. We will in the following refer to "heavy mesons" as " B" (where we mean the 
realistic example of a B c meson, or the toy model for a B q meson) and "light mesons" as "7r" (where 
the realistic example is r] c , and the toy- model application would be the pion or also the kaon for 
mi / m 2 ). 

2.2 Definition of LCDAs for light pseudoscalar mesons 

Following [1,2] we define the 2-particle LCDAs of a light pseudoscalar meson via 



(7r(P)|gi(y) [y,x]~f^ 5 q 2 (x 
(<K(P)\qi(y) [y,x] i^q 2 {x 
(n( p )\Qi(y) [y,x]a /ll/ j 5 q 2 (x 



-if* 



du e i ( u P-v+ u P- x ) 



mz 



Vti 4>*(u) + — — Zf, g n (u) 
2p • z 



du e ^P-y+up-z) ^(u), 



du e i i u P-y+ u P- x ) ^ a ^ u ) 



2D 



(3) 



with two light-like vectors z^ = y^ — x^ and p^ = P^ — m 2 /(2P ■ z) z^, and u = 1 — u denoting the 
light-cone momentum fraction of the quark q\. The gauge link factor is denoted as 



[y, x\=V exp 



1 9: 



[ dt{y-x)- A(ty + (1 - t)x) 
Jo 



(4) 



4>-k{u) is the twist-2 LCDA, while 4> p (u) and (p a {u) are twist-3. For completeness, we have also 
quoted the twist-4 LCDA g n (u) which, like the 3-particle LCDAs, will not be considered further in 



this workU All LCDAs are normalized to 1, such that the prefactors in ([3|) are defined in the local 
limit x — > y. In the definition of cp a (u), we have included a factor 2>/(D — 1), such that the relation 
between fj, n and p, n from the equations of motion (see below) is maintained in D / 4 dimensions. 

2.2.1 Equations of motion 

The equations of motion (eom) provide relations between the matrix elements defined in ©. Fol- 
lowing [2] we obtain 
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D - 1 



In — 1 

(2 - D) 4> a (u) H — 4>' a {u) 



4>n(u) + g n (u) = (mi + m 2 ) 4> P {u) + ■ ■ ■ 
(mi + m 2 ) 4>-k(u) + • . • , 



(2u - 1) (^p(it) + ^ - 4>' a (u) = (mi - m 2 ) <M«) + • ■ • , (5) 



where the ellipsis denote contributions from 3-particle LCDAs which we do not specify here. In the 
local limit the contributions from the 3-particle LCDAs drop out and integration of ([5]) yields 
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mi 



and 



Hit = , Att = Mvr ~ (™1 + m 2 ) (6) 

mi + rri2 



u (f)p(u) = — H . (7) 



o 



2 2 M7r 



Notice that the relations (|5I6I7[) hold for the bare (unrenormalized) parameters and distribution 
amplitudes. 

2.2.2 Tree-level result 

At tree level, and in leading order of the expansion in the relative velocities, the quark and the 
antiquark in the NR wave function simply share the momentum of the meson according to their 
masses, pf ~ mj/(mi + m2)P^. For "light" mesons this implies^ 

</) w (u) ~ <j) p (u) ~ g n {u) ~ S(u - uq) , (8) 

with uo = mi/(mi +m2) and uo = m 2 /(mi-|-m 2 ). Consequently, all positive and negative moments of 
the distribution amplitudes are simply given in terms of the corresponding power of uq. In particular, 
the Gegenbauer coefficients are given by 

an = + 3 ? ; f du Mu) (2u - 1) 2^n+^ ^ 

3(2 + n)(l + n) y 3(2 + n)(l + n) 

Notice that /2„- ~ at tree- level and the corresponding LCDA <f> a (u) can only be determined by 
considering the corresponding one-loop expressions (see below). The tree- level solutions © fulfill 
the eom-constraints from ([S]). 



2 Notice that there are additional two-gluon LCDA for flavour singlet mesons which we will not consider here, 
because in the non-relativistic limit glueballs decouple from the qq states and the 2-gluon LCDA is only generated by 
higher-order relativistic corrections. For the definition of the 3-particle LCDAs, see [1,2]. 

3 This behaviour can also be obtained from the "dense medium limit" in the instanton model [30]. 
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2.3 Definition of LCDAs for heavy pseudoscalar mesons 

We define the 2-particle LCDAs of a heavy pseudoscalar B meson following [3,4] 

if B (ji)M 



(0\(q)p(z) [z,0] (h v ) a (0)\B(Mv)) 



2 



75 



, (io) 



J a/3 



where is the heavy meson's velocity, t = v ■ z and z 2 = 0. Here fg is the (renormalization-scale 
dependent) decay constant in HQET. The Fourier-transformed expressions, which usually appear in 
factorization formulas, are given through 



dw e~ iujt (j)±(oj) , 



(11) 



where u denotes the light-cone energy of the light quark in the B meson rest frame. 



2.3.1 Equations of motion 

The equations of motion again provide relations between different LCDAs. Including the effect of 
the 3-particle LCDAs 4^4, Xa, Ya as defined in [31] (see also [32,33]), we derive 

D — 2 f u 

u (j) B {u) - m 4>%{u) + — ^— / dr] [^(77) - ^(rj)] 

ru; poo Jt Q 

= (D-2) / dr, / -f- [*a(v,Z)-*v(t,,Z)] ■ (12) 
Jo Juj-i) K °K 

The relation (|12p is trivially fulfilled at tree-level and we will show below that it also holds after 
including the a s corrections to the NR limit. In [31], Kawamura et al. discuss a second relation 
which in the massive case reads 

(uj + m) <j>g{(^) + (u! — 2 A — m) (f>^(uj) 

1 n, r f ^ aM+XaM] -, (d - 2] f in r f m 

duJ Jo Ju>- V ? JO Ju-r, £ d£ 

with A = Mb — rn^. We will show below that the equation (|13j) does not hold beyond tree level, since 
the integral on the right-hand side involving our result for the 3-particle LCDA Xa does not converge. 
This confirms the criticism raised in [24, 34] that (j 13f) is not consistent, since the renormalization 
prescription of light-cone operators in HQET and the expansion into local operators do not commute. 
Notice that in contrast to (|12p . the derivation of (|13p involves derivatives with respect to z 2 7^ 0. 

If one neglects the 3-particle distribution amplitudes in (|12l) , one arrives at the so-called Wandzura- 
Wilczek relation which has first been discussed for a massless light quark in [4]. The generalization 
to the massive case reads 

/•w 2 
J o dr] [4> B {ri) - 0+(r?)] ~ — - [w^(«) - m^+(w)] , (14) 

which again holds for the bare parameters and LCDAs in D 7^ 4 dimensions. 



2.3.2 Tree-level result 

By the same arguments as for light mesons, at tree- level the quark and the antiquark in a heavy 
meson just share the total momentum according to their masses, such that uj = m. In the NR limit, 
the 2-particle LCDAs of a "heavy" meson are thus given by 

<f%(w) ~ <f>g(w) ~ 6(u - m) . (15) 
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Moreover, at tree level, the moments of the heavy meson LCDAs can be related to matrix elements of 
local operators in HQET [3]. The zeroth moment (0|#75 h v \B) = —ifsM determines the tree-level 
normalization of the distribution amplitudes 4>%{t = 0) ~ 1. For the first moment, one has the 
general decomposition 

(0\(q)piE»(h v ) a \B(v))^--^- [{av^ + brf) (1 + ^)75]^ ■ (16) 

Multiplying by (757 A1 ) ( g Q and taking into account the finite light quark mass in the NR set-up, the 
equation of motion for the light quark implies a + 4b = m. The equation of motion for the heavy 
quark is obtained by multiplying with v 1 *, from which one obtains a + b = A, independent of the 
light-quark mass. This implies 

4A — m , A — to 



From this we can read off the first moments at tree-level 

% * — 4A — Tft 

(u)+ ~ T — (0|g 7 5 4- (in-D) h v \B) = a = — — - , (17) 
f B M 3 

i *~ 9 A -I- m 

(0\q l5 r/L + (in-D)h v \B) = 2b + a = —^, (18) 

f B M o 

where we introduced the light-like vectors = z^/t and n+ = 2v^ — n_. In the non-relativistic 
limit A = m, and we obtain 

(uj)± ^ m . 

Notice that the light-quark mass drops out in the sum 

+ = 2A. 

We stress that the relation between moments of 4>^(u>) and local matrix elements in HQET does not 
hold beyond the tree-level approximation [23,24,35]. 



3 Relativistic corrections at one-loop 

The NR bound states are described by parton configurations with fixed momenta. Relativistic gluon 
exchange as in Figure [2] leads to modifications: First, there is a correction from matching QCD 
(or, in the case of heavy mesons, the corresponding low-energy effective theory HQET) on the NR 
theory. Secondly, there is the usual evolution under the change of the renormalization scale [5,23]. 
In particular, the support region for the parton momenta is extended to < u < 1 for light mesons 
and < u < 00 for heavy mesons. In this section we collect the results for LCDAs for "light" and 
"heavy" mesons including the first-order matching corrections from relativistic gluon exchange 

<PM = ^ + ^^)+0(al). (19) 

3.1 Light mesons 

3.1.1 Local matrix elements 

We first consider the leading-order relativistic corrections to the local matrix elements which are 
given by the vertex-correction and the wave-function renormalization of the quark fields. We find 



fir fir 



NR 



47r V mi + m,2 m<i 



(20) 
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Figure 2: Relativistic corrections to the light-cone distribution amplitudes. The dashed line indi- 
cates the Wilson line in the definition of the LCD As. 
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. (21) 



where m n ~ m° s + to™ in the on-shell scheme. Our result for the decay constant is in agreement 
with [36] and the results for fi n and p, w are consistent with the eom-constraints in (|6|). 

3.1.2 The twist-2 LCDA 4> n (u) 

Let us start with the case of equal quark masses, e.g. in case of a non-relativistic r] c bound state, 
which may also serve as a toy-model for the pion LCDA0 

The first-order relativistic corrections arise from the collinear gluon exchange diagrams in Figure O 
where we also have to take into account the wave-function renormalization of the external quark lines 
(see Appendix [A] for details). The local limit of the light-cone matrix element ([3]) determines the 
relativistic corrections to the NR decay constant (|2(jp (in this case, the diagrams with the gluon 
attached to the Wilson-line do not contribute). The remaining contributions to the NLO correction 
for the leading-twist LCDA contain an UV-divergent piece, 



#(* 



I div. 



duV(u,v)<f><®(v), 



which involves the well-known Brodsky-Lepage evolution kernel [5] 

Mm, , / 1 



V(u,v) 



1 + 



v 



u 



U -6(v 

V 



The finite terms after MS-subtraction read 

,2 



In 



+ 4 



7712(1/2 

u(l — u) 



u 



u + 1 + 



1 + 



V 



(22) 



(23) 



1/2 



u 



u 9(1/2 -u) + (u++u) 



(1/2 -«) 

Here the plus-distributions are defined as 



( du { . . . } + f(u) = J\u {...} (f(u) - /(1/2) 
f du { . . . } ++ /(«) = J du {...} (f(u) - /(1/2) - /'(1/2) (u - 1/2) 



(24) 

(25) 
(26) 



4 We should keep in mind, however, that typically non-perturbative analyses from lattice QCD and sum rules find 
pion distribution amplitudes that are broader than the asymptotic one, while the non-relativistic model assumes LCDAs 
which are narrower. Therefore the application of the toy model to the very pion case should not be taken too seriously. 
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Table 1: The moments (£ n )7r at the non-relativistic scale ji = m. 



n 


2 4 6 8 10 


NR limit 

NLO 423) (for a s = 0.2) 
4r © (for < R = 0.2) 



0.067 0.011 0.004 0.002 0.001 
0.067 0.008 0.001 0.000 0.000 



From this it follows that 



du (p^ (u; fi) = I duu <ffl (u; fi) 



such that the general normalization conditions f du<f> n {u) = 1 and f duu(f> v (u) = 1/2 are not 
changed. Furthermore, our result for the distribution amplitude obeys the evolution equation 

^(u- fi) = / dv V(u, v) ^{v- fj) + 0(a 2 s ) . (27) 



dh\fl 7T 







An independent calculation of the leading- twist LCDAs for the r/ c and J /ip meson has been presented 
in [25]. Our result is not in complete agreement with these findings. In particular, we find that the 
LCDA quoted in [25] is not normalized to unity as it should be. 

On the other hand, at the non-relativistic scale fi ~ m, the distribution amplitude shows a 
singular behaviour at u ~ uq = 1/2. As a consequence, the convergence of the Gegenbauer expansion 
is not very good at the non-relativistic scale, with the Gegenbauer coefficients a n in ^ only falling 
off as l/\/n (and alternating signs). A better characterization of the LCDA at u ~ m is given in 
terms of the moments 

(DM = I du (2u - l) n ^(n; fi) , (28) 
J o 

which are linear combinations of Gegenbauer coefficients of order < n. This corresponds to an 
expansion of the LCDA in terms of a delta-function and its derivatives, 

M«i/i) = 2X;^-i (B) (2 U -l) (DM- (29) 

n=0 

Results for the first few moments {£ n ) n are shown in Table Q] for the strict non-relativistic limit, 
including the NLO corrections from (|24p and comparing with the non-relativistic corrections of order 
discussed by Braguta et al. in [26]. Keeping first-order corrections in only, this formally 
amounts to the replacement 

^ R (u) = 5{u - 1/2) + ^ 5"{u - 1/2) + 0(v* R ) . (30) 



In particular, this fixes the moment (^ 2 )- K = The authors [26] propose a resummed formula, 



/NR/. a 1 n( 1 - U NR \ n I 1 + ^NR 



The comparison in Table \T\ shows that for ~ a s (m) ~ 0.2, the effect of the corrections is 
qualitatively and quantitatively similar to the a s corrections from (J2^ 
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It is also interesting to determine the correction to the first inverse moment of the LCDA which 
appears in QCD factorization formulas 



u 



3 2.73 + 1.08 In 



in- 



Finally, we quote the result for the derivative of 4>ir{u) at the endpoints 



0Uo;^) = -0Ui;m) 



4tt 



12 In 



A* 



in. 



+ 0(a 



s I - 



(32) 



(33) 



which is sometimes discussed in the context of non-factorizable contributions to hard exclusive reac- 
tions [37,38]. 

For non-equal quark masses, the NLO corrections to the MS-renormalized twist-2 LCDA are 
given by 



+ 4 



In 

™ ' K 

( u(l — u 



m 2 (u Q - u) 2 



1 



1 + 



1 



Uq — U / Uq 



It 



6(u - u) + 



\ (uq - u) 2 

The first moment now becomes 

o. s Cf 



+ 2 5'(u - u ) ( 2n (l - n ) In 



»0 



1-Uq 



n <— > u 

Uq <-)• U 

- 2n - 1 



dn n (it; u) = uq + 



4tt 



4 

q lll_ 2 2 

3 «n ml- 



7(1 - n ) 



In Ur 



38 



Uq - [U <-> n j 



(34) 



(35) 



3.1.3 2-particle LCDAs of twist-3 

The twist-3 LCDAs for the 2-particle Fock states are obtained in the same way as the twist-2 one. 
After absorbing the corrections to the local matrix elements into the renormalized values for fj, n and 
fin, we obtain a UV-divergent piece 



f'Mldi 



div. 



1 + 



1 



6(u - n) + 



n <-> n 
n <-> n 



no — n x 

and a finite NLO contribution to the twist-3 LCDA associated to the pseudoscalar current 



(36) 



In 



+ 4n (l - n ) 



1 



1 



1 + 



1 



Uq — U 



9(u - u) + 



n <-> n 
n <->• n 



(n - n) 2 

In particular, the first moment of (f> p (u) now reads 



+ <5'(n — no) In 



no 



1 - n 



f 2n - 1 
1 ("o - «) 



o 



(in n </> p (n; u) = no H — — 

4tt 



3 In + 6no In n — 4 ) no — (n <-> no) 



m 



A" 



which is in agreement with the eom-constraint from (|7|). At the endpoints we now have 



</> p (0;n) 



a s C F ( 2 + 2n 
47T \ no 



/' 



In ^ -2 +C>K) 



(37) 



(38) 



(39) 



and similar for </> p (l;u) with mi <-> m.2, i.e. no «-> no. For the twist-3 LCDA associated to the 
pseudotensor current (whose normalization factor starts at order a s ), we simply have 



<t> a {u) 



u 



(n - n) + 



n <-> n 
n <->• n 



+ e>K). 



(40) 



In contrast to the other 2-particle LCDAs in 0, we find that <^ CT (n) is not given by a delta-like 
distribution in the NR limit and has support for < n < 1. 
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3.2 Heavy mesons 

The calculation of the LCDAs for a B c meson (which again can be considered as a toy model for 
LCDAs of B q mesons with m^ m„) goes along the same lines as for the r\ c case. However, important 
differences arise because the heavy 6-quark is to be treated in HQET which modifies the divergence 
structure of the loop integrals (notice that in our set-up, a charm quark in a B c meson is treated as 
"light"). As a consequence, the evolution equations for the LCDAs of heavy mesons [23] differ from 
those of light mesons. 



3.2.1 The LCDA <£+(. 



oj 



Let us first focus on the distribution amplitude 4^(oj) which enters the QCD factorization formulas 
for exclusive heavy-to- light decays. In the local limit we derive the corrections from soft gluon 
exchange to the decay constant in HQET. We find 



NR 



/m(m) = 



1 + 



Att 



31n-^-4) +0{aj) 



m 



(41) 



Notice that the decay constant of a heavy meson exhibits the well-known scale dependence in HQET 
[39]. The remaining NLO corrections to the distribution amplitude <P~b(oj) contain an UV-divergent 
piece (details of the derivation can be found in Appendix |B|) 



0b ( w ;a0 
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2oj 
e 



6(m — oj) 0(oj — m) 
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8(oj — m) 



1 1 



1 - In 



m- 



(42) 



and a finite piece 
4 + ' 1} (u;;m) _ 



OJ 



In 



(oj — m) 1 



(m — oj) 9(oj — m) 



m(m — oj) oj{oj — m) 



+ 4 



8(2m - oj) 
(oj — m) 2 



++ 



+ 



4 6(oj- 2m) _ 5(oj - m) (\ 2 
(oj — m) 2 m \2 m? 



hi 



H 2 3tt 2 



+ 



+ 2 



(43) 



with an analogous definition of plus-distributions as in (|25)26p . Notice that, in order to separate the 
UV divergence coming from the longitudinal momentum integration, we have introduced an auxiliary 
parameter [if = 2m to split the support region of the LCDA into two parts. The distribution 
amplitude in (|4^|) obeys the evolution equation 



d 



din fi 



4>%(uj;ij) 



a s CF 
Ait 



doj'^\u,u';fj,) <p+(oj' ]f i) + 0(a 2 s ) 



(44) 



where the anomalous dimension ^\oj,oj'; fi) can be read off the UV-divergent terms in (|42p and is 
given by [23] 



7 i 1} (oj, oj'; /x) = (r« sp In £ - 2)8(oj - oj') - T% oj 



oj 



6(oj' -oj) 9(oj - oj') 



oj'(oj' — oj) oj(oj — oj') 



(45) 



with riulp = 4. 



In contrast to the light meson case, the normalization of the heavy meson distribution amplitude 
is ill-defined. Imposing a hard cutoff Atjv 3> m an d expanding to first order in m/A\jy, we find 



doj (fi~^(oj; fi) ~ 1 



Ait 



2 A 2 r A? 



UV 



IT 2 ' 

+ 12 



+ 0(a 2 s ) + 0(m/Au V ) 



(46) 
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and similarly for the first moment 

Auv 



1 



Auv 



duo lu <f)^(u>; n) ~ 



a s CF 
4tt 



2 In 



o 



A 2 



+ 6 



+ 0{a 2 ) + 0(m/A 



uvj 



(47) 



The last two expressions provide model-independent properties of the distribution amplitude which 
have been studied within the operator product expansion in [35]. Our results are in agreement with 
these general findings. 

We finally quote our result for two phenomenologically relevant moments in the factorization 
approach to heavy-to-light decays [24, 35] 







0t(w;/i) 



— (l- ° scf 

m \ 4-7T 



l,,|i 2 . /U 2 3vr 2 „' 
- In 2 - In + — - 2 
2 m z m z 4 



+ 0(a2), (48) 



and 



a B {ji) = agfa) = In ^ + ^ [8(3] + 0(a 2 ) , 
m 47T 

where (j = is the Riemann zeta function and we defined 



(49) 



OJ 



OJ 



(50) 



The leading-order scale-dependence of these quantities is in general given by 



dxl_ 

dlnfi 
ddB 
dlnfi 



a,Cp 



1 + 



An 



(51) 
(52) 



In particular, in the non-relativistic limit cr^ (//) = (<7B(/z)) n and therefore the a s correction on the 
right hand side of (JMJ) does not depend explicitly on In//. For arbitrary values of n, we find for the 
scale dependence of the logarithmic moments 



da 
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dlnfi 



na 



(n-l) + a sC F p(i) 



4n 



cusp 
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(n+1) 



[n/2] 



+ 2n! 
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J (n-2j)! 



+ 0(a 



2 ) 



(53) 



where [x] denotes the greatest integer less than or equal to x. 
3.2.2 The LCDA <j) B {u) 

A similar analysis can be performed for the other 2-particle LCDA of the B meson. We now obtain 
for the UV-divergent piece (see also Appendix |B|) 
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(54) 
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The finite contributions read 



4> B ' 1 \u; f i)=2 

+ 2uj 
+ 4m 



In 



A* 



(u — m) 2 
2 



In 



(u; — m) 2 
w — 2m) 
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a; (a; — m) 
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+ 4m 



6»(2m - uj) 
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m z 



(55) 



The distribution amplitude in (|55p obeys the evolution equation 
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din 



(j> B (uj;(i) 



o 



4"7T 
47T 

(1) 



do;' 7i 1) (w,(j / ;^) <j) B {u)'\n) 



(56) 



where the anomalous dimension kernels j 1 ^' (uj,uj'; fi) and 7_i(u;,u/; /i) can be read off the UV 
divergent terms in (|54p (see Appendix [B] for details) 



7i 1) (w,w'; / u) = 7^ i; (u>,a/;A t ) — F 



,(i) 0(u)'-u) 



cusp 



UJ' 



cusp 
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(57) 
(58) 



Among others, the knowledge of 7_ is essential to check the factorization of certain correlation 
functions appearing in sum-rule calculations for B — > it form factors within SCET [40]. 

Another new result are the first positive moments of the LCDA (f> B (co) as a function of a hard 
cutoff Atjv ^> m i 



Ait 



2 A 2 ,, A 2 ,, T 12 
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Auv 



"■UV 
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+ 0(a 2 ) + 0(m/A uv ), (59) 



Auv 
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A 2 



+ 2 



+ 0(a 2 ) + CJ(m/A 



(60) 



which are again expected to be model- independent. Actually, these moments can already be derived 
in the Wandzura-Wilczek approximation. From the solution of (|14p in D = 4 — 2e dimensions 



'B 



(w; /i) = (1 - e) / dr? 



m 



we obtain the bare (unrenormalized) moments 

A»+V5(A) + (l-e) (w»>+ 



n + 1 — e 



(61) 



(62) 



which result in the same MS-subtracted moments as in (|59l60p . i.e. the 3-particle LCDAs only 
contribute subleading terms to these moments in our case. 
We finally quote the quantity 



^-( 0;/ u) = ^lln4 + O(«s), 
B Ait m m 2 8 

which plays a role in sum-rule calculations for heavy-to-light form factors [40,41]. 



(63) 
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3.2.3 3-particle LCDAs and equations of motion 



In order to verify whether the equations of motion (|12|13p hold after including first order relativistic 
corrections, we have to compute the 3-particle LCDAs which arise at order a s in the non-relativistic 
limit. Without going into details, we quote our results for the bare LCDAs that enter (|12I13|) 



cx s Cf 8(u> — m + £) 
47T 2m 

a s Cj? 5(u — m + £) 
47r 2m 

o s cf r / 1 

47T \ £ 



~ + ln 72 
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(4m - 30 



£ C 2 



2m 



1 , M 
- + ln-^ 1 

e 77V 



+ 1 m d 5(e 



(64) 



We show in Appendix O that the eom-constraint (|12p is indeed fulfilled to order a s . On the other 



hand we find that (|13p does noi hold beyond tree level since the ^-integral involving our result for 
the 3-particle LCDA Xa is ill-defined for £ — > oo. Since we again expect the radiative tail of the 
3-particle LCDAs (which determines the large-£ behaviour) to be model- independent, the failure of 
([TBI beyond tree level should be considered a general feature. 



4 Renormalization-group evolution 

In physical applications, the light-cone distribution amplitudes are required at the hard-scattering 
scale /x which is set by the momentum transfer in the exclusive reaction. The evolution from the 
"soft" scale m to the hard-scattering scale fi resums large logarithms In fi 2 /m 2 . In this section we 
study the evolution of the NR distribution amplitudes to leading logarithmic (LL) approximation. 
For "light" mesons we focus on the twist- 2 LCDA <fi n (u) and for "heavy" mesons we consider the 
2-particle LCDAs 4>^(uj) and <j)g(uj). 

4.1 The twist-2 LCDA (f> w (u) 

The evolution of the twist-2 LCDA cj) n (u; /i) is described by the Brodsky-Lepage kernel (|23p . To solve 
the evolution equation (|2"T)) . one projects the distribution amplitude onto Gegenbauer polynomials 
which are the eigenfunctions of the evolution kernel, 

Mu; H) = Gnu ^1 + a n {n) C z J 2 (2u - 1) j . (65) 

The respective coefficients are obtained from ([9]) and have the LL evolution 

/ n \ -7n/A> / 9 Tl+1 1 \ 

\<Xs(lM))J \ (n + l)(n + 2) j^j J 

with /?o = (33 — 2n/)/3 (for illustration, we will use = 3 in the numerical examples). For very 
large n>lwe have 7„ ~ —4 Cp In n , which implies that the effect of higher Gegenbauer coefficients 
becomes less important at high scales. 

We show in Table [2] the LL evolution of the moments (£ n ),r defined in (128p . starting from the 
tree-level result in the NR limit, ^(u; /j,q = m) = 5(u — 1/2), for two values of 77 = a s (n)/a s (m) and 
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Table 2: The moments (£ n ) n as a function of the evolution parameter 77 = a s ()i)/a s (m). 



n 


2 4 6 8 10 


LL (r, = 1/5) 
LL (77 = 1/25) 
asymptotic 


0.126 0.048 0.025 0.015 0.010 
0.173 0.070 0.038 0.024 0.016 
0.200 0.086 0.048 0.030 0.021 



in the asymptotic limit. In contrast to the moments {£, n ) n , the phenomenologically important 1/u 
moment is a linear combination of an infinite number of Gegenbauer moments, 



00 



(u- 1 ) n (p)=3}2a 2:i (fi). (67) 

j=0 

In order to study the evolution effects from the non-relativistic scale, where (u~ 1 ) n ([io = m) = 2, 
towards fj, 3> m, it will therefore be crucial to control the effects of higher Gegenbauer coefficients. 
For this purpose, we find it convenient to consider model parameterizations which are obtained from 
a slight modification of the strategy developed in [42]. Our ansatz involves three real parameters 
a > 0, b > 0, < t c < 1, 

Mu) = r( " - ; / dti-lnty- 1 f(2u-l,it 1 / b ) + f(2u-l,-it 1 / b ) (68) 
r(a;-lnt c ) J I J 

with r(a; b) = J fe °° dt t a-1 e -i and the generating function of the Gegenbauer polynomials, 

. 00 



Cl' 2 {2u-1), (70) 



Performing the t- integration in (J68J), one finds 



4>-k{u) = 6im 



00 

cos 



(t) r( " , ^-^,""' ) ^ + " 



n=0 

from which one reads off the Gegenbauer coefficients a dd = and 

(_l)«/2 r(o; -(l + n/6) lnt c ) 



for n even. (71) 



{n/b + l) a T(a;-\nt c ) 

For t c — > our ansatz reduces to the asymptotic distribution amplitude and for t c — > 1 it is equivalent 
to one of the models discussed in [42], where the Gegenbauer coefficients show a simple power-like 
fall-off (in this case with alternating signs). As observed in [42], for values of a < 3, the model induces 
some pathological behaviour at u = 1/2. In our ansatz this is regularized by the cut-off parameter 
t c < 1. The qualitative behaviour of the Gegenbauer coefficients for large n now depends on t c : 

• For moderately large values of n, we have 

l«n«n crit E-6(l + l/lnt c ) : \a n \ ~ (n/6)~ a , (72) 
i.e. a power-like fall-off with n as in [42]. 

• For asymptotically large values of n, one obtains 

6(-lnt c )° _1 t n J h 

ra>n cri t : \a n \ ~ — ^—^ — , (73) 

1 (a; — mi c ) n 

i.e. an exponential fall-off with n, which renders the contribution of very high Gegenbauer 
coefficients irrelevant. 
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We now reconsider the evolution of the tree- level result in the NR limit ^(u; liq = m) = 5(u— 1/2) and 
fix the model parameters (a, b, t c ) in (|68p from the first three non- vanishing Gegenbauer coefficients 
using ([!]), 

a ~ 3962 

a 2 (m) ~ -0.5833 ^ fe „ Q ^ 

a 4 (m) ~ +0.4583 V f „ Q g9g3 (74) 

a 6 (m) ~ -0.3906 I 

V ' ' n crit H49 

The fact that a < 1 and n cr it 3> 1 reflects the bad convergence of the Gegenbauer expansion in the 
NR limit. Still, the model parameterization reproduces the Gegenbauer coefficients with n <^ n cr j t 
and the value of the first inverse moment (u -1 )^ to a very good accuracy (see the first line in Table[3|). 

The same strategy can be applied for scales fx > m. The LL evolution towards larger scales 
depends on r\i = a s ((Xi)/a s (m). For illustration, we consider r/i = 1/5 and rj2 = 1/25 and obtain 

a ~ 1 2679 

a 2 ( m ) ~ -0.2160 | b „ a87Q8 

04(^1) - +0.1079 J> ^ ^ ^ 98n (75) 



a 6 (/ii) ^ -0.0679 
and 



"-crit. ^ 45 



a ~ 2.1451 
-2 ^2 - -0.0800 | ^ 

a 4 (^ 2 ) ^ +0.0254 J- =>■ ^ „ Q 941g (76) 



oeCAfc) - -0.0118 



"-crit. ^ 14 



We observe that the parameter a increases under evolution, which is related to the growth of the 
anomalous dimensions for larger values of n, leading to a steeper fall-off of the Gegenbauer coefficients 
at larger scales. Effectively, for moderately large values of n, one has 

a(m) « a(m) - Inn, . (77) 

The parameter b is only slightly increasing while t c is decreasing under evolution. The critical value 
of n is quickly decreasing from n cr it — 1149 at fx = m to n cr it — 14 at u = u 2 . Figure [3] shows the 
evolution of the model LCDA as a function of u. For r/ = 1/5 the functional form still "remembers" 
the non-relativistic profile, while for 77 = 1/25 it is already close to the asymptotic form. Table [3] 
compares the first few Gegenbauer coefficients using the exact projection of the delta-function and 
the model parameterization (|68p . We see that the differences are tiny and the model gives a good 
approximation. We also quote the result for the first inverse moment which slowly evolves from the 
NR value, {u~ 1 ) n (fXo = m) = 2, towards the asymptotic value (u~ 1 ) 7r (fx — ► 00) = 3. We clearly 
see that the model gives a better description for relatively low scales than a truncated conformal 
expansion (|65p with n < 6 (i.e. with the same number of input parameters as our model). The latter 
can be improved, however, by considering the averaged moment 

^ / 3iriax 3 max 1 \ 

(w^m - 2 H a M + Yj a M > ( 78 ) 

\j=o 3=0 J 

which accounts for the alternating sign behaviour of the Gegenbauer coefficients. Using this improved 
truncated conformal expansion, we find that the 1/u moment is given by (n- 1 )^ = (2.04,2.57,2.82) 
at respective scales (m, fx±, 1x2), which is very similar to the predictions of the model parameterization 
(see Table EJ). 
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Table 3: LL evolution of the first few Gegenbauer coefficients starting from the NR distribution 
amplitude 4 > tt{ u ) — S(u— 1/2). For each value of the evolution parameter 77 = a s (/i)/a s (m) we show 
the results for the exact expression, the model parameterization (|68p and a truncated conformal 
expansion (|65[) with n < 6. We also quote the value for the first inverse moment (u -1 )^. 





a.2 


0,4 


o 6 


a 8 


aio 


012 


Oi4 


016 




exact (77 = 1) 
model 

conformal exp. 


-0.583 


0.458 

* 


-0.391 

* 


0.346 
0.346 


-0.314 
-0.314 
(trun 


0.290 
0.289 
cation n 


-0.271 
-0.269 

< 6) 


0.255 
0.253 


2.00 
1.45 


exact (77 = 1/5) 
model 

conformal exp. 


-0.216 


0.108 

* 


-0.068 

* 


0.048 
0.048 


-0.036 
-0.036 
(trun 


0.029 
0.028 
cation n 


-0.023 
-0.023 

< 6) 


0.020 
0.019 


2.55 
2.47 


exact(?7 = 1/25) 
model 

conformal exp. 


-0.080 


0.025 

* 


-0.012 

* 


0.007 
0.007 


-0.004 
-0.004 
(trun 


0.003 
0.003 

cation n 


-0.002 
-0.002 

< 6) 


0.002 
0.001 


2.81 
2.80 
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Figure 3: Approximation of the NR distribution amplitude 4>ir( u ) = H u ~ 1/2) in terms of the 
model parameterization (thick solid line) and its evolution for 77 = 1/5 (dashed line) and 
77 = 1 /25 (dotted line) . The asymptotic LCDA is shown for comparison (thin solid line) . 



4.2 The LCDA <j>+(u) 

The evolution of the LCDA (/^(ui; /i) for scales m < /j, < M is described by the Lange-Neubert kernel 
(|45p Fl The solution of the evolution equation (|44p can be written in closed form as [35] 



T{2-g) f°°du/ , 



W > Y W < TP fl O O U< 

— — 2-Fi 1 - 9, 2 - g; 2; — 

fJLQ J LOy \ u> 



where uj < = min(w,u/) and lo > = max(w,w'). The evolution is controlled by the functions 



V = V(fi,lMi) 



with r cusp ~ aft, 7 



2tt ' ^ ~ 



da 



■Tcusp (oi) 



da' 



/9(a 



n + 7(a) 



2 ^ (we use n/ = 4 in the numerical examples) and 

« s (m) 



r cU s P (a) 2Cf , a s (//o) 
g = g(jjt, po) = I da „ v ■ ~ — — In ■ 



a s (A»o) 



/3(a) 



A) a a (/x) 



(79) 



(80) 



(81) 



^Notice that above the &-quark mass scale one should match the LCDAs in HQET onto LCDAs in QCD, see e.g. [43] 
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Figure 4: Evolution of the heavy meson LCDA ^(uj; fi) starting from the tree-level result 
Ho = m) = 6(u) — m), where we assumed a s (m) = 1 (thick solid line). The curves (dashed, 
dotted, thin solid line) correspond to r\ = a s ([i,)/a s (m) = 1/2, 1/5, 1/10, respectively. 




The hypergeometric function 2-^1(0, b; c; z) has the series expansion 

^ T(a + n)T{b + n)T{c) z n 

2 F 1 (a,b;c;z) = > vFTvTF? T S f 

^— ' r(a)r(6)r(c + n) n! 

n=0 

Starting from the tree-level result in the non-relativistic limit (/)^(co; /j,q = m) = 5(lo — in), we obtain 
for scales [i > m the relatively simple expression 

,v-2 7ES £(2-5) /w>y W< 9 _„.o.^< 



m</>+(u,; M ) = e ^7*ff v*i ^ _^ ^ 1 - 5 , 2 - 5 ; 2; , (82) 

tree 1 (g) V 171 J <jJ > V U > / 

where now = min(cj,m), o; > = max(w, m), (7 = g(fi,m) and V = V(fJ,,m). Fixing the value of 
a s (m) at the NR input scale, we may study how the shape of (j)^(uj; fj) is changed by evolution effects. 
In Figure H] we have plotted (|52"j) for a s (m) = 1 and three different values of rj = a s (fj,)/a s (m). As 
expected, the evolution drives the initial delta-function shape towards a flatter distribution. In the 
double- logarithmic plot on the right-hand side in Figure IH we may read off the asymptotic behaviour 
of 4>b{io\ fi) for u — > and u> — > 00. As argued on general grounds [23], the LCDA develops a linear 
behaviour for uj — > 0, whereas for a; — ► 00 it tends to fall off slower than 1/u at higher scales. This 
can also be seen by comparison with Figure where we plot the evolution of another LCDA with 
initial condition (f>^(to; hq = m) = iv/m 2 e~ u} ^ m . 

In Figure [6] we show the corresponding evolution of the phenomenologically relevant moments 
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Ab(m)/A b (/x ) 

4 - 
3 
2 





0.10 0.15 0.20 0.30 0.50 0.70 1.00 




■AO 0.15 0.20 0.30 



0.50 0.70 1.00 



Figure 6: The moments Xb{^)/ Ab(mo) & n d as {(J-) — ln////xo as a function of the evolution parameter 
rj = a s (fJ,)/a s (fJ>o)- The solid line refers to the initial condition (j>~^(u>; (j,q = ni) = 8{lo — m) and the 



dashed line to <j}%{u>; jJLo 



ui/m 



j / m 



where we assumed a s (m) = 1. 



As (n) and ctb(h) denned in (|48I50|) . From ([79j) we find the closed formulas 

i r(i-g) r*" (<»\ 9 ^ ^ 

3(/*) r (!+5) Jo w \/V 



Ab(a») 

<7 B (/i) = g(l-g) 4*3(1, l, l - 5, 2 - 5; 2, 2, 2; 1) 

F 2 (l- (7,1-^,1-^2,2-^1) 



1-5 



W£(w;/z ) / — — ln-^+(a;;/io). (84) 



w V/W J Jo ^ V/W A* 

In general, the evolution of the moments A# and erg thus depends on the shape of the LCDA 
4>^{uj) [23]. For our examples, ^(u; /j,q) = 5{uj — fJ-o), respectively ^(u; (Mq) = uj/^q e~ w l^ , the 
a; integration can be performed explicitly, leading to relatively simple analytic expressions. It is 
also possible to approximate the factors (u/ fio) 9 = 1 + g ln(u;///o) + • • • In this approximation, the 
evolution for the moment As(/x) can be entirely determined in terms of Xb{i^o) and gb^o)-, see 
also [24]. 

4.3 The LCDA (f) B (u) 

The evolution of the LCDA <j)^{oj; fJ,) is somewhat more involved, because of the possible mixing 
with the 3-particle LCDAs. In addition, for a non-vanishing light quark mass m / 0, we have seen 
in ([56]) that the LCDA <p^(uj;fi) mixes into <pg(u;; fi). In the following, we concentrate on possible 
applications in realistic B q decays (where we can set m = 0), neglecting the contributions from 3- 
particle LCDAs which is left for future work. In this approximation, the solution of the evolution 
equation 

- 0b( w >'M) = _ - / duJ ' (t>s(u;';n) + 0(a 2 s ) 



dlnfj, ay ^' Air J 



^ J terms with m/0 (85) 
1 contributions from 3-particle LCDAs 

can be obtained in a similar way as for ^(lo; fi) [23,35]. The details of the derivation can be found 
in Appendix [Dl As a result, the solution for ^(to; fj,) can be written as 

- e^'S^-s) r—KMlM>) (-Y 2^(1-9,1-9,1,^). (86) 

r (5) Jo ^> \wy v 
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m4> B {u,n) 
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to/m 
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Figure 7: Evolution of the LCDA <fig(to; n). Notations and conventions as in Figures HISI The initial 
conditions at the scale /io = m are 4>g(to] fJ-o) = 5{to — m) (upper row), <fig(to; no) = 9{m — to)/m 
(middle row) and </>^(w;/io) = l/me~ w / m (lower row). 



In Figure [7] we illustrate the evolution of 4> b {uj;^) for three different initial conditions at the scale 
fi = m: 

• 0b(w;A*o) = 5{uj - m), 

• (f) B (u;a ) = 0(m-to)/m, 

• cj>-{u;u Q ) = I/me""/" 1 . 

The first example corresponds to the strict non-relativistic limit (where the neglect of the light quark 
mass in the evolution equation may be considered as inconsistent). The second and third example 
follow from the Wandzura-Wilczek relation (|14p for the initial LCDAs (j)~^{to] considered in the 
previous subsection. While the behaviour of (j) B (to;fi) at small values of to depends on the model 
for the initial distribution, the radiative tail for large values of to is again universal. More precisely, 
the solution ([86]) of the (approximate) evolution equation suggests that 4> B (to; fi) also falls off slower 
than 1/to at higher scales, while the slope of the LCDA at to = tends to vanish under evolution, 
independent of the initial behaviour of the distribution amplitude. 
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5 Summary 



Non-relativistic qq bound states have been used as a starting point to construct light-cone distribution 
amplitudes for light mesons in QCD and heavy mesons in HQET. At the non-relativistic scale, the 
leading 2-particle distribution amplitudes can be approximated by delta functions, fixing the light- 
cone momenta of the quarks according to their masses. After including radiative corrections from 
relativistic gluon exchange, the distribution amplitudes cover the whole physically allowed support 
region, < u < 1 for light mesons and < uo < oo for heavy mesons. In this paper, explicit 
expressions for 2-particle distribution amplitudes of twist-2 and twist-3 for "light" mesons (with 
quark masses m\ ~ 771,2) have been calculated to first order in the strong coupling constant. In the 
same way, next-to-leading order expressions for the 2- and 3-particle distribution amplitudes have 
been derived for "heavy" mesons (where m\ S> 777,2). We also studied the evolution of the 2-particle 
distribution amplitudes under change of renormalization scale. 

Our results apply to the physical situation of a hard exclusive reaction, that involves bound states 
of heavy bottom or charm quarks, with large momentum transfer, for instance, B c — > r] c £v [27,44-47], 
e + e~ — > J/iprjc [48-51] or 7*7 — > r) c [52]. Moreover, from the divergence structure of our explicit next- 
to-leading order results, we could derive certain model-independent properties which also hold for 
bound states of relativistic quarks. In this way we used our calculation as a toy model to derive new 
results for the B meson distribution amplitude 4>g, as the cut-off dependence of positive moments, 
the anomalous dimension kernel and the solution of the evolution equation in the Wandzura-Wilczek 
approximation. The toy model also allowed us to address an issue that has been controversial in the 
literature, i.e. the question if the constraints from the equations of motion hold beyond tree level in 
the heavy meson case. 
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A One- loop corrections to (f>n(u) 

We briefly summarize our results for the individual diagrams in Figure [2] in the light meson case (in 
Feynman gauge). For simplicity we present the results for the leading- twist LCDA (j) n (u) and stick 
to the case mi = 777,2 — m - 



A.l Vertex diagram 



Starting from the NR limit 4>tt(u) = 5{u — 1/2) and performing the loop-integral in D 
dimensions, one obtains for the first diagram in Figure [2] the distribution 



I a (u) oc 4r(e) 



m 2 (1 - 2u) 2 



4u 2 - 4u - 1 
(1 - 2u) 2 



[u6(l - 2u) + u6(2u - 1)] 



2c 



(87) 



The integral contains an UV-divergence reflected by T(e). The IR-divergence at u = 1/2 can be 
isolated with the help of a plus-distribution which we introduce via (|26p . With this we obtain 



I a {u) oc 4 



— h In - r 
e 777/(1 — 2u) z 



4u 2 - Au - 1 
(1 - 2u) 2 



[u9(l - 2u) +u6(2u- 1)] 



( -+3 1n^-2)<5(n-l/2). 
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Notice that the term with 5'{u — 1/2) vanishes due to the symmetry u <-> u in the equal mass case 
(the "++" -distribution actually coincides with the usual "-(-"-distribution in this case). The local 
term determines a correction to the decay constant and does not contribute to ^(u). 



A. 2 Wilson-line diagrams 

For the second diagram in Figure [2] we obtain 



rl/2 

I b {u) oc 8T(e) / dv 
Jo 



m 2 (l - 2v) 2 J 2v-l 
Convoluting with a regular test function, we get 



[5(u - 1/2) - 5(u-v)]. 



i ,i 
duf{u)Ib{u) oc / duf(u) 











8it0(l-2u)/l n 2 
l-2u \e + m 2 (l-2u) 2 



(89) 



(90) 



The other Wilson-line diagram in Figure [2] is obtained from I b by symmetrization 



u — > u. 



B One- loop corrections to 4>b(uj) 

In the following we present our results for the diagrams in Figure [2] in the heavy meson case (in 
Feynman gauge). We compute the first order corrections to the NR limit at the matching scale 
fx ~ m starting from 4> B (u)i n ) = 5(uJi n — m), and the general anomalous dimension kernels related 
to the renormalization of 4>~^{uj; fi) and 0^(o;;/x). The latter are extracted from the UV-divergent 
parts of the diagrams, where we consider arbitrary input functions <p^{uj- m ) and also keep track of 
the light quark mass m ^ 0. Notice, that a possible mixing of the 3-particle LCDAs into (fig(u);n) 
is not considered. Some care has to be taken when performing the collinear limit (which amounts to 
setting the transverse momentum of the incoming light antiquark to zero). 



B.l Vertex diagram 

The loop-integral in the first diagram of Figure [2] reads 

I a( U; )\ fj fun <5(w - Win + n_J) 



^ oc - J duj in J [dl] 



I-{u))J J m J 1 J [vl + i0][l 2 + i0][{l- k) 2 + i0] 



n_l-u in + -m fc 



k±4j „ ; m? i , . / k±-l 



rn p n + L - — + u- 




(91) 



where uj- m = n^k with being the momentum of the incoming spectator quark and k^ — P is the 
spectator-quark momentum after the interaction with the gluon. We also performed the collinear 
limit k± — > 0, which requires to keep terms of order k± ■ l_\_/k\. 

Let us first consider the fixed-order corrections to the NR limit, where <^(ui n ) = 5{uj- m — m). 
Then the loop integrals simplify according to 

oc [\dl]- r . no 2m ~Z t1 , 5(u-m + n-l). (92) 

aV ; J 1 1 [vl + iO][l 2 + iO][P-2mvl + iO] y ' y ' 
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Performing the loop-integrals in D = 4 — 2e dimensions, one is left with the distributions 

2 6(m — uj) 9(uj — m) 



I+(u) oc 2wr(l + e) 



(m — uj) 2 J \(m — ui) 2 m(m — uj) uj(u) — m) 



Auj 



9(2m - uj) 
(m — uj) 2 



+- 



A 6(oj - 2m) n 

+ 4uj - 2u 

(m — uj) z 



9(m — uj) 6(uj — m) 



m(m — uj) lo(u> — m) 



( - + hxt--4) 8(u-m), 
e m z 



(93) 



I a (uj) oc 



2(1 - e) r(e) / 



m 

+ 2mr(l + e) 



(m — uj) 2 



6(m — uj) 



9(m — uj) 6(uj — m) 



(m — uj) 2 J I (m — uj) 2 m(m — uj) m(ui — m) . 





+ In 


\ ^ 1 


(m - 


-uj) 2 _ 


- 2 


9(m 
m - 


-co)' 

- UJ 


+ 



1 



(m — uj) 



+ 4m 



2uj 



m 

6(uj — m) 
uj(uj — m) 



0(2m - uj) 
(m — uj) 2 

,2 
^2 



+ 4m 



+- 



6(ui - 2m) 
(m — uj) 2 



+ 2 ( - + In ^1 - 4 ) S(u - m) 
1 e m z 



(94) 



The integration over uj determines the local vertex correction to be absorbed into the decay constant 

(95) 



Focusing now on the UV-divergent contributions to the integration kernels in the general case, 
we find 



9(co in - uj) 



<j>z(uj- m ) + O(e ) 



(96) 
(97) 



Notice that only the kernel in I~ receives an UV-divergent piece, which can be traced back to the 
appearance of a factor (n+l) in the numerator of (|91|) . The fact that the kernel of l£ is UV-finite is 
in line with the findings of [23]. 

B.2 Wilson-line coupling to heavy quark 

In this case there is no mixing between fig and (f>~^ since the light-quark propagator is not involved 



I b (uj) oc - / duj- m / [dl] 



5(uj — Ui n + n-l) — 5(uj — uji n ) 
(n-l) [vl + iO] [I 2 + iO] 



0|(Win) • 



(98) 



Inserting the non-relativistic LCDAs and performing the (n + l) and l± integrations, one is left with 
the parameter integral (k = — n,-l) 



I±(uj) oc 2T(e) 



dk 



fi 2 e 7E \ e 5(uj — m — k) — 5(ui — m) 



k 2 



k 



(99) 



Notice that the remaining integral induces an additional UV-divergence, which has to be isolated by 
introducing appropriate plus-distributions. We find 



I b (to) oc 2uj 



- + In 

e 



(uj — m) r " 



6(uj — m) 
to(u — m) 



1 1 u 2 1 2 M 3tt 2 \ _ f 

-~ + - m — k + - In — H — o(uj - m). 

e 2 e m 2 2 m 2 4 / v 

(100) 
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The UV-divergent contribution for l£ corresponds to the result for the diagram (Dl) in [23] (with 
cj' = m). The UV-divergence from k — ► oo is a peculiarity of the heavy meson wave function. It 
is related to the cusp-anomalous dimension involving the heavy quark (characterized by a time-like 
vector v^) and the soft Wilson line (characterized by a light-like vector n^_). The resulting 1/e 2 terms 



are universal for and (j) B , 



div. 



X 



11, P 

-~ + - m — 



duj\ 



9(UJ - W in ) 



U3 — OJh 



(101) 



B.3 Wilson-line coupling to light quark 

In this case the loop integrals mix cj)g into 4>g (but not vice versa), 

5(u> — u>i n + n_Z) — 8{u — oj\ a ) 



oc 2 / duj 



dl) 



ri-l 



(n^l)[(l-k) 2 +i0] [l 2 + i0] 







m ■ 



k i -l i 



fc2 W in ^ fc2 

Inserting the non-relativistic LCDAs and performing the (n+Z) and Z^ integrations, we find (A; = n_Z) 

m — k ( fi 2 e lE \ e 5(k — m + to) — d(ui — m) 



1 



^(^in) 



(102) 



I+(lo) oc 2T(e) / rf/c 
'o 



k 2 



k 



2lo 



+ ln 



A* 



(w — m) 5 



0(m — u) 
m (m — to) 



A* 



+ | - + 21n^ + 4 U(w-m), (103) 



7~(w) oc 2T(e) / dA; 

'o 



jjL 2 e 1E \ e 6(k — m + to) — 5(u> — m) 



1 



+ In 



A- 2 

2 



(u> — m)' 



8(m 



m — oo 



The UV-divergent contributions to the integration kernels are identified as 

u 9{oj\ n — uj) 



ltd 



2 2 
''Idiv. K 7^bH + - 



cZcji, 



2 2 



Win (Win - W) 



+ - / duJ in 



m 8(uj 



Win (Wi, 



[<^(w in )-0+(u,)]+0(e°), 

[«/>i(Win) - 

in) _ (/) + M ] +O (6 ). 



(104) 



(105) 



(106) 



Our result for Ij is in line with [23]. In particular, there are no additional UV-divergences related 
to cusp anomalous dimensions since the light-quark and the soft Wilson line are characterized by 
the same light-cone vector n/\ For a non- vanishing light-quark mass, the result for I~ implies that 
the LCDA 4>~^ mixes into cftg under evolution. In the massless case, however 
independently (at least to leading logarithmic approximation). 



, an d 4>b evolve 



C Equations of motion for heavy meson LCDAs 

In this appendix we show that the eom-constraint (I12p holds after including first order relativistic 
corrections to the NR limit. We first evaluate the right-hand side of (|12p using our explicit results 
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for the 3-particle LCDAs from (f64"|) 

(D-2) d V -f -[y A ( v ,0-*v(v,0} 



47T 



1 Ll Z 

— h In — k ) m 5(u — m) 



2uj 1 



+ to- 
rn V £ ym — ujy 



+ 1 +21n 



(m — a;)" 



6(m 



For the expressions on the left-hand side of (|12|) . we obtain 



(m os - m barc ) S(u - m) + 



_barcN , a sCF f 2cJ ( 1 ( /T 



+ In- 



( - + 21n ^ 

(cj — my 



An \m \e (m — uj) 2 

2 
2 



1 0(m-u) 



2 ) 6(u - m) - ( - + 2 In At + 4 ) m 5(w - m) 
1 e 777/ 



and 



D-2 



dr) [<f>%(rj) - <P B {r])] 



4tt 



m \£ (m — LV) z 



+ 2 In 



im — uj)' 



9(m — uj) 



+ I - + 21n 

— mj z 



2 }6(u - m) 



Noticing that 



r m os _ m b are) = ^Cf (^ + 3ln E_^ 



2 ) 

s 1 1 



(107) 



(108) 



(109) 



47T \e 

we see that the equation of motion (fT2|) is indeed fulfilled after including the a s corrections. 



D Solution of RGE for 4> b (uj) in WW approximation 

We derive the solution of the evolution equation (|85[) for the S-meson LCDA (J) B (u>;fi), ignoring 
the possible mixing with 3-particle LCDAs and neglecting the light quark mass m. We follow the 
analysis in [35], where a closed form for the LCDA <p B (uj; fi) to LL approximation has been given. 
When the 3-particle LCDAs are neglected, 4> b (uj; /i) can be related to ^(uj; fi) by the Wandzura- 
Wilczek relation (|14p . This is also reflected in the leading-order result for the anomalous dimension 
kernels j_ (uj, uj'; /i) from (I57j) and j^\uj, uj'; /x) from (|45j) . Noticing that 



-UJ 



r« p in £-2) 8( u -n), 



;uo) 



and 



/>oo 

/ d?7 f(rj) 
Jo 



d 






f 


duj 


Jo 


-f 


d.r 


Jo 


X 



*» dJ 



6{uj - uj 1 

UJ — UJ 1 

-x))-f{u)\ 



UJ 



dy 



d f 1 dx 
duj J x J 1 - 



/Ki-y)) 



dr] f(rj) 



UJ — Tj 



111) 
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and 



dv f(v) j-^T" f ~ 
o I du Jo f] 



uj 9(u' — uj) 
uj' {uj' — to) 



UJ 



00 dx 
x(l + x) 



;/( w (i + x))-/H] = / d v /( v ) 



d r°° dx f x dy 
duj Jo x(l + x) Jo 1 + y 

to 9{r] — uj) 



fWi + v)) 



rj (77 - u) 



we find that the anomalous dimensions fulfill the relation 



d p dw (i) , (i) 
duj Jo V 



;ii2) 



(113) 



Therefore, the functions $> b (uj; /i) = uj d<j) b (uj; n) / du) and 4> b (uj; fj) obey the same evolution equation 
to LL approximation. Using an intermediate result from [35], we write the solution for & b (uj;[i) as 



Mo/ 



JU m=l v 



(-l) m+1 r(l + m- g) 



m + <?)r(l + m)r(m) 



/ </; \ — m / </; \ rn—q 

x \me{uj-uj') y—^j +{g-m)e{J-u) 



(114) 



with V and g from (|8QI8ip and we assumed < g < 1. We finally perform the u;-integral to obtain 
4> B (uj;[i) from $ b (uj; /j,), and the summation over m which leads to hyper geometric functions with 
the result, 



mo7 



\ 9 



x ^(w-w') - 2Fx(l-g,l-g,l ) u'/uj) 



UJ 



;ii5) 
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